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1. INTRODUCTION 

Supercapacitors (SCs), also called ultracapacitors or electric double layer capacitors, are electrical 
energy storage technologies. Compared to batteries, they have a better specific power density, a longer cycle 
life (over 104 cycles) and a better efficiency [1]. They are also characterized by lower series resistance and a 
more significant equivalent capacitance [2], [3]. These strengths make them attractive for several applications 
such as electric or hybrid vehicles for peak power [4]-[6], wind power generation [7], solar power generation 
[7], [8], or simply consumer electronics [9]-[12]. They can be used independently, or in combination with 
other energy storage technology such as lead—acid batteries [13], lithium batteries [5], [14] or hydrogen fuel 
cells [15], [16]. 

During their operation, supercapacitors are subject to aging, which causes a progressive degradation 
of their performances. There exist some factors that can accelerate aging mechanisms such as high voltage, 
high operating temperature, and voltage imbalances that take place when several supercapacitors are 
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connected in series [17]. The industrial expansion of supercapacitors requires a deep understanding of aging 
mechanisms and the development of techniques for the estimation of their state-of-health (SoH) during 
operation. Such techniques are very important to predict and anticipate supercapacitor breakdown. 

During SC operation, the online identification of its internal resistance and capacitance is important. 
SC degradation is characterized by an increase of the internal resistance or a decrease of the equivalent 
capacitance [17], [18]. Therefore, these two parameters constitute good indicators of SC state-of-health. 
There identification can facilitate the planification of preventive maintenance tasks and improve system 
safety. 

SC resistance and capacitance can be determined based on the impedance measurements using 
electrochemical impedance spectroscopy (EIS) technique [19], [20]. However, EIS is an offline 
measurements technique, which means that the system’s operation must be interrupted. The objective of the 
present paper is to develop an online (real-time) identification technique. 

Gualous et al. [21] proposed an identification algorithm based on extended Kalman filter (EKF) to 
estimate buffered energy of a SC used for a solar application. The SC was modeled by three-branch circuit. 
They assumed that RC circuit parameters are constants with the aging time. In [22] also proposed an online 
identification technique based on the EKF to estimate the temperature and the state of charge of a 
supercapacitor. In [23] proposed two real-time observers for the prediction of supercapacitors’ capacitance 
and the resistance. One is based on EKF and the other is based on interconnected observers. Both observers 
used RC circuit to model supercapacitor. They assumed that capacitance and resistance remain constant. The 
results obtained by the two online observers are then compared with an offline characterization. 

Nadeau et al. [24] developed a technique for the real-time identification of SC resistance and 
capacitance for vehicular applications. In order to model SC, they chose RC circuit in which the capacitance 
is assumed to vary linearly with bias voltage. Sliding mode observer with an adaptation gain was used to 
identify SC internal resistance and capacitance. Based on a Lyapunov-based adaptation law, developed a 
real-time resistance/capacitance identification method for supercapacitors [25]. Compared to other methods 
in literature, the stability and the convergence of their method is guaranteed by Lyapunov’s direct method. In 
[26] proposed an online identification scheme based on particle filter (PF) for the estimation of state-of- 
charge, SoH and temperature of SC by combining electrical equivalent circuit model and thermal model. 
They found that the proposed estimation scheme performs SoH estimation with acceptable accuracy. In [27] 
propounded a nonlinear state-space model for the online estimation of capacitance and resistance of SC. The 
proposed model takes the capacitance variation and self-discharge effects into account. The estimation error 
achieved with the proposed model is approximately 5%. According to [28] suggested a procedure integrating 
fuzzy logic and artificial neural network (ANN) to estimate supercapacitors’ internal resistance and 
capacitance. ANN can be a powerful technique to provide a robust identification for systems that are 
subjected to uncertainties [29]. However, this technique usually requires long offline tests to have the enough 
training data. 

The EKF becomes a recognized and a widely used technique for state estimation in several 
nonlinear dynamic systems. One of the advantages of EKF in our case is that it only requires terminal 
measurement of voltage and current, which can be done during SC operation. However, the main difficulty 
using EKF in our case consists in setting initial values of the SC parameters. These values must be close to 
real values. Otherwise, the EKF is unable to accurately identify the correct values of SC parameters (state 
vector) [21]-[23]. This limitation was also highlighted in this paper (section 5). Unfortunately, it is not 
always easy to accurately estimate the initial system state. In order to address this limitation, a 
complementary PID corrector was integrated to EKF scheme in the present paper. In literature, some authors 
have already investigated the combination between EKF scheme and PID corrector in order to improve 
systems control [30]-[32]. This combination allows better stability and better convergence [30]-[32]. The 
novelty of the present paper is to investigate the use of EKF combined with PID corrector to identify SC 
parameters. 

Section 2 presents the experimental procedure used to characterize SC. In section 3, the equivalent 
circuit model used in our method is presented. In section 4, SC state vector is determined based of the chosen 
equivalent circuit model. In section 5 the extended Kalman filter scheme is presented. The complementary 
PID corrector is then detailed in section 6. In section 7, the obtained results are reported and discussed. 
Finally, conclusion and some future works are recommended in section 8. 


2. EXPERIMENTAL STUDY 

In this experimental study, two supercapacitors maxwell technologies BCAP1500 and BCAP350 are 
subjected to charging and discharging cycles under the two current profiles shown in Figure 1. The applied 
currents correspond to the real currents during SC operation. A schematic of the experimental bench is 
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presented in Figure 2. Both bias voltages response (Figure 3) and applied currents are measured during 
experiment and saved using the acquisition board NI cDAQ 9178 and LabView software. Sampling time 
used for voltage/current measurements is set to 0.01s. The experimental bench used is presented in Figure 4. 
In the sections that follow, the proposed identification method is applied to the two supercapacitors. 


Current profile applied to SC BCAP1500 Current profile applied to SC BCAP350 
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Figure 1. Current profiles applied to the two SC, (a) BCAP1500; (b) BCAP350 
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Figure 2. Schematic of experimental bench used for SC characterization 
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Figure 3. Voltage response profiles of the two SC, (a) BCAP1500; (b) BCAP350 


The internal resistance and capacitance of a SC are affected by operating temperature as 
demonstrated in [24]. It is thus important to conduct the experiment in a constant temperature. For this 
reason, the two SCs were placed inside a temperature-controlled chamber in which the temperature is set to a 
constant value of 65° C. 
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Figure 4. Experimental bench used for SC characterization 


The experimental results obtained at this stage will be used to identify SC parameters for the two 
SCs using the new online identification method proposed in this paper. It is worth noticing that 
voltage/current measurements can be done in real time and does not require the interruption of SC operation. 

In the second part of this experimental study, the same experimental bench described in the present 
section is used to apply Maxwell characterization test to the two SCs. Maxwell test is an offline identification 
test that allows the determination of internal resistance and capacitance. Although maxwell test provides 
accurate and reliable results [33], it requires the interruption of system operation and therefore, cannot be 
performed in real time. In the present study, the values of resistance and capacitance obtained by maxwell 
test are used as reference values to compare the results obtained by the proposed real-time identification 
method and to assess the accuracy of this method. 
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Figure 5. Maxwell identification test 


The maxwell characterization test is based on the assumption that the SC can be modeled by a 
simple RC equivalent circuit. More details on this model are given in section 3. Based on this model, the 
voltage U, measured on the terminals of the SC can be described as the sum of the voltage on the terminal of 
a resistor and a capacitor, denoted respectively U, and U}. The SC is submitted to the current profile shown 
in Figure 5. According to Maxwell test procedure [30], and based on voltage response (Figure 5), internal 
resistance and capacitance are determined using (1) and (2): 


R=2 (1) 
d 
_ Iata 
C= (2) 


The different parameters used in (1) and (2) are illustrated in Figure 5. The internal resistances 
obtained using maxwell characterization test are respectively R = 6.67 x 1074 Ohm and R = 3.32 x 
107’ Ohm for supercapacitors BCAP1500 and BCAP350. The capacitances obtained using Maxwell 
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characterization test are C = 1498.21 F for supercapacitor BCAP1500 and C = 349.56 F for the 
supercapacitor BCAP350. 


3. SUPERCAPACITOR MODELING 

In order to model the electrical behavior of a SC, several equivalent electrical circuit models have 
been proposed in literature. The most common are simple RC model [24], [34], Zubieta model [35], simple 
pore model [10], CPE model [36], multi-pore model [37], and fractional-Order Model [38]. Increasing circuit 
sophistication always leads to a more complex formulation and generally requires high computational 
resources, making it more difficult to identify SC parameters in real time [39], [40]. Zubieta model [35], 
presented in Figure 6, is a simple model that evaluates the overall real behavior of the supercapacitor. It is 
composed of a series combination of a resistor R and a capacitor which is characterized by variable capacity 
over time C. The capacitor represents the canonical capacitance effect of SCs. The resistor basically 
represents the electrolyte and electrodes resistances. 


i(t) R oO) 
U, (t) U.(t) 
U,(t) 


Figure 6. Equivalent electrical circuit model 


i(t) is the charge/discharge current. The bias voltage U,(t) on the terminals of SC can be described 
in (3): 


U,(t) = U,(t) + U,(t) = R.i(t) + eae i(t). dt (3) 


It has been proven that the capacitance C(t) varies linearly with the voltage U,(t) [35], [41], [42]. 
Therefore, the capacitance C(t) is expressed as a function of voltage U,(t) using the following linear 
expression: 


C(t) = Co + a.U,(t) (4) 


Co is the capacitance value when the voltage U,(t) is null and a is the proportional coefficient. In the 
rest of this paper, Cy and a are assumed to be constant. 


4. EXPRESSION OF SC STATE VECTOR 
Extended Kalman filter (EKF) allows state estimation of nonlinear dynamic systems during each 
time period. In our study, the state vector chosen in our work is given in (5). 


Get eal bas he 5 
l=? (5) 
X,) la 


EKF is used to estimate, at each time period, the state variables of a continuous nonlinear system 
defined by (6): 


ie = A.X(t) + w(t) = f(X, i) (6) 


Y(t) =G.X(t) + v(t) = g(X ©), i) 
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Y(t) is the system output, which corresponds to the bias voltage U,(t) in our study. t is the time 
index, X(t) is the state vector, w(t) is the state noise and v(t) is the measurement noise. The determination of 
w(t) and v(t) are detailed in section 5. In this section, the determination of A and G is detailed. 

We assume that, during a cycle of several charges/discharges, the variation in resistance is 
negligible. We also consider that the capacitance increases with an almost constant rising slope as a function 
of the voltage, with a proportional coefficient a. We can thus assume that the variations of resistance R as 
well as of the proportionality factor a are zero as a function of time: 


dR 
aa (7) 
da 
ao (8) 


The derivative of the voltage U, at the terminals of capacitor is described as shown in (9): 


auc es a.Uc—-C . 
dt (a2.u2—c?) 9) 


Based on the expression given in (4) and (9), capacitance derivate C can be described as shown in (10): 


dc a auc a= a.(a.U¢ — C) . 
dt ` dt (a?.U¢? — C?) ` 


(10) 


Finally, based on previous equations, the derivate of state vector can be expressed as shown in (11): 


X = A.X (11) 
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Where: 


—i(t) 
X4".X14°-X3" 
X,.i(t) 
X4".X14°-X3" 
Xe I(t) 
X4".X4°-X3" 
—X,.i(t) 
X4 X1 -X3 


(12) 


The estimated system output, which corresponds to the voltage at the terminals of the SC, is denoted U, (t). It 
can be expressed as a function of state vector X as shown in (13): 


U,(t) = R.i + U, = G.X = [1i 0 0].X (13) 


5. EXTENDED KALMAN FILTER SCHEME 

One of the advantages of using extended Kalman filter (EKF) is its ability to deal with nonlinear 
systems, as it is the case capacitance and the bias voltage behavior of SCs. It has been used with success in 
order to estimate systems states in many applications [16], [18], [19], [29]. Its recursive structure and low 
computational cost make it suitable for real-time identification of system state parameters. EKF allows the 
estimation of system state, represented by the state vector X in our case, at each time step (in real-time). In 
addition, EKF only requires terminal measurement of voltage and current, which can be done in real-time. A 
functional diagram of EKF is presented in Figure 7. 

The EKF is used to estimate the state vector X of a system described by (6). In this equation, Y(t) is 
the system output (voltage at the terminals of the SC). X(t) is the state vector presented in (5). t is the time 
variable. The time increment used for our simulation is set to 0.01s. That corresponds to sampling time used 
for experimental voltage/current measurements. w(t) and v(t) represent respectively state and measurement 
noise. In our study, it is assumed that state and measurement noises are white, gaussian, uncorrelated and 
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centered with zero expectation. It is also assumed that the noises w(t) and v(t) are known respectively by 
covariance matrices W and V. The accurate determination of W and V terms is difficult since these noises are 
related to unmodeled dynamics and uncontrolled variability in parameters. In [24] conducted an empirical 
study using offline datasets to estimate the variance terms, which allows the estimation the terms of 
covariance matrix W and V. The same procedure is used in the present study to determine W and V. 
Therefore, the covariance matrix of noise condition is as shown in (14) and (15): 


W = diag(0.17 0.00017 0.00012 0.12) (14) 
V = 0.01 (15) 


The EKF is based on two important steps: prediction and estimation. The first step consists on predicting the 
state vector X and the covariance matrix of the estimation error P(t), which satisfies as shown in (16): 


AT 
P(t) =E |(x - 8). (x-2) | (16) 
Where, X is the predicted state vector. From (16), a Riccati differential equation is derived to compute P(t): 
P(t) = F.P(t) + P(t). FT — P(t).GT.V7t.G.P(@ +W (17) 


The transition and observation matrices used in the extended kalman filter are defined as the following 
Jacobians (the first order partial derivative of vector functions): 


_ af 

E (18) 
_ ag 

G= ak (19) 


The second step of the extended Kalman observer is the correction step, where the Kalman filter 
gain K is calculated based on the covariance matrix of the estimation error P(t). Using this gain, the predicted 
state vector X(t) and the covariance matrix of the predicted measurement error P(t) are updated, which leads 
to the estimated state X(t) and P(t). Figure 7 shows the functional diagram of the extended Kalman filter. 
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Figure 7. Functional diagram of extended Kalman filter 


Before using the proposed EKF scheme, it is important to verify the observability of the system. 
Because of the non-linearity of our system, a linearization around an operation point is required [43], [44]. 
The observability is then verified around an equilibrium point X, using the Jacobian matrix J(X,) as defined 
in [45]. After the calculation of J (Xe), we found that rank (J (X,)) = 4, which correspond to the order of our 
system. Therefore, our system is observable around its equilibrium point Xe [45]. 
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In the rest of this section, the results obtained using simulation are presented, the limitation of using 
EKF is explained to highlight the need for complementary correction. then, the proposed PID corrector is 
presented in the next section. First, simulations have been conducted using conventional EKF scheme 
explained in the present section (without any correction) for the two supercapacitors and using the cyclic 
current profiles described in Figure 3. In order to estimate the resistance and capacitance using conventional 
EKF algorithm, experimental bias voltage is also used (Figure 4) with the same time increment (0.01s). The 
obtained resistance and capacitance values using conventional EKF algorithm were compared with those 
obtained using Maxwell test, considered as the reference values. Figure 8 and Figure 9 show the obtained 
results. As shown in these figures, simulations have been performed using different initial values of 
capacitance and resistance. 


Resistance profile for SC BCAP1500 Resistance profile for SC BCAP350 


Simu. 1: initial resistance = 0 mO Simu. 1: initial resistance = 0 mQ 
Simu. 2: initial resistance = 1 mQ Simu. 2: initial resistance = 1 mO 
Simu. 3: initial resistance = 5 mQ Simu. 3: initial resistance = 5 mQ 
Resistance of Maxwell test = 0.667 mQ Resistance of Maxwell test = 3.32 mQ 
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Figure 8. Estimation of supercapacitor’s resistance using different initial values for two SC, 
(a) BCAP1500; (b) BCAP350 


Figure 8 (a) and Figure 8 (b) shows that, regardless of the initial value chosen, the resistance 
converges toward values close to experimental resistance (Maxwell test). The absolute errors between 
simulated resistance and Maxwell test resistance are respectively 3.71% and 4.02% for BCAP1500 and 
BCAP350. However, as shown in Figure 9 (a) and Figure 9 (b), the capacitance profile is highly dependent 
on the initial value. When this later is far from the real value of capacitance, EKF is unable to converge to 
this real value. As the capacitance is one of the key parameters to diagnosis supercapacitors state-of-health 
[17], [18], it is very important to accurately identify this parameter during SC functioning. In the rest of this 
section, the EKF scheme is analyzed to understand why the obtained capacitance is far from the correct 
value. 
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Figure 9. Estimation of supercapacitor’s capacitance using different initial values for two SC, (a) BCAP1500; 
(b) BCAP350 
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As explained in section 5, the EKF determines SC state parameters by two important steps: 
prediction and correction. During prediction step, state parameters are predicted recursively as shown in 
equation (6). At each time increment t,, the derivate X, of state parameters is first calculated using equation 
(6). Then, using Euler's method as shown in (20), the state vector X, at time step t, is predicted using the 
derivate X, and the state vector X,., determined at time step t,.;: 


LAA FACA (20) 


Where At is the time span. 

Since state vector is calculated recursively, the predicted state values at each increment are very 
dependent on the initial values. Therefore, the slightest error in initial state vector Xp) will be propagated to 
the other state vectors X,,...,X,. Unfortunately, it is not always easy to accurately estimate the initial system 
state. Most of the time, the correction step proposed in EKF scheme (section 5) allows to correct the 
prediction errors associated to misestimation of initial state vector. This is one of the benefits of correction 
step, proposed in EKF scheme. Using the gain K, the state vector Xp predicted at the first step of EKF 
scheme can be corrected as shown in (21). This includes the prediction errors associated to misestimation of 
initial state vector. 


Xo = p +K x Err (21) 


Where Xc is the corrected state vector and Err is the error between simulated bias voltage U, and 
experimental bias voltage Y. The error Err is determined as shown in (22): 


Err = (U, — Y) (22) 
In our case, the simulated bias voltage U, is obtained by: 


U. 
U, = G.X = [1 i 0 0]. : =U, +R.i (23) 


a 


As can be seen in (23), among the four state parameters U., R, C and a, only supercapacitor voltage U, and 
equivalent resistance R are used to estimate bias voltage U,. The other parameters C and a are not used. 
Consequently, the correction term K X Err depends only on U, and R. Once the two parameters U, and R are 
corrected, the error Err tends towards zero. Therefore, the correction term K x Err also tends towards zero. 
As a consequence, the parameters C and a are no longer corrected. This explains why any error in initial 
value is propagated to all the other values. 


6. COMPLEMENTARY PID CORRECTOR 

For an effective capacitance correction, the bias voltage simulated during correction step must be 
determined considering the value of capacitance. Based on Zubieta model (Figure 6), the voltage on the 
terminals of SC can be expressed by: 


Usc(t) = J, ~i(t).dt + Uo + R.T (24) 


Where Up is the voltage on the terminals of the capacitor when t = Os. Considering time discretization, the 
bias voltage can be approximated by: 


1. : 
Uscltn) = Lizo,...n g ilti). At + Uo + Rn- itn) (25) 


Unlike the bias voltage expressed in (23), as shown in (25) takes into account the capacitance. Therefore, the 
bias voltage expressed in (25) will be used in the complementary correction of capacitance. 

A proportional—integral—derivative corrector (PID corrector) is a control loop that use feedback 
results to correct system output. Because of its high efficiency, PID is a widely used corrector in industrial 
area. It is very appropriate to system requiring continuous correction, which corresponds to the present case. 
At each time increment, the error value, expressed in (26), is calculated as the difference between simulated 
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bias voltage U,,(t,,) and experimental bias voltage Y(t,). It is worth noticing that bias voltage U,,(t,) is 
estimated in this case using (25). Therefore, the capacitance is considered in the estimation of bias voltage. 
Based on this error, a correction is applied to the capacitance at each time increment based on proportional 
(P), integral (I) and derivate (D) as expressed in (27). 


e(tn) = U, (ty) = Y (tn) (26) 


e(t;)—e(t;_1) 
Corrected = Coredicted + Kpe(tn) + Kı diro,...n €(ti)-At + Kp i At — (27) 
Cyredicted represents the capacitance obtained in prediction step of EKF. Coorrectea represents the 
corrected capacitance. Kp, K; and Kp denote the coefficients for the proportional, integral, and derivative 
correction terms respectively. Figure 10 represents the functional diagram of extended Kalman filter 
combined to PID corrector. 
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Figure 10. Block diagram of extended Kalman filter combined to PID corrector 


7. RESULTS AND DISCUSSION 

Using the new algorithm described in Figure 10 (EKF scheme combined to the complementary PID) 
and the cyclic current/voltage profiles presented above, simulations have been executed for the two different 
supercapacitors BCAP1500 and BCAP350. As the first simulations, different initial values of 
capacitance/resistance are tested. The results are presented and discussed in the present section. In Figure 11 
(a) and Figure 11 (b), a comparison is made between the simulated bias voltage obtained by the new 
algorithm and the experimental bias voltage. We can clearly see that the two profiles remain very close 
throughout the simulation for the two SC. The maximum errors between simulated and experimental bias 
voltage for BCAP1500 and BCAP350 are respectively 0.007V (0.3%) and 0.01V (0.4%). These results 
clearly demonstrate the ability of the new algorithm to accurately simulate the bias voltage. 

The real-time evolution of capacitance obtained using the new algorithm is depicted in Figure 12. In 
this figure, the capacitance obtained by simulation is compared to the capacitance obtained using Maxwell 
test, considered as the reference value. As can be seen in this figure, the simulation was performed using 
different initial values of capacitance in order to assess the ability of the new algorithm to deal with errors 
associated to misestimation of initial capacitance. From Figure 12 (a) and Figure 12 (b), it can be clearly seen 
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that, regardless of the initial value used, the capacitance converges towards the value obtained by Maxwell 
test, contrary to the same simulations carried out without PID corrector (Figure 9). This convergence is the 
result of the PID corrector integrated to the EKF scheme. At the end of simulation, the errors between 
capacitance obtained by the proposed method and obtained by Maxwell test are respectively 1.02% and 
2.87% for BCAP1500 and BCAP350. 
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Figure 11. Experimental and estimated bias voltage profiles for two SC, (a) BCAP1500; (b) BCAP350 





Capacitance profile for SC BCAP1500 Capacitance profile for SC BCAP350 
Simu. 1: initial capacitance = 100F 1750 Simu. 1: initial capacitance = 100F 
Simu. 2: initial capacitance = 1000F Simu. 2: initial capacitance = 230F 
Simu. 3: initial capacitance = 1500F 1500 Simu. 3: initial capacitance = 350F 
Simu. 4: initial capacitance = 2000F Simu. 4: initial capacitance = 460F 
Capacitance of Maxwell test = 1498.21F 71250 Capacitance of Maxwell test = 349.56F 
D 
ia 000 
g 
a 750 
oO 
O 
500 
250 J 
0 
0 50 100 150 200 250 300 350 0 250 500 750 1000 1250 1500 1750 2000 
Time (s) Time (s) 
(a) (b) 


Figure 12. Estimation of supercapacitor’s capacitance after using PID corrector for two SC, 
(a) BCAP1500; (b) BCAP350 


By comparing the two Figures 9 and 12 the same initial values are used for both simulations. The 
capacitance results obtained using EKF show a divergence from the Maxwell value which depends perfectly 
on the initial value proposed however, by analyzing the results of the simulation using EKF with PID 
corrector; we see that the capacitance converges towards the optimal value (Maxwell capacitance) and does 
not depend on the initial value. We note from results obtained in sections 5 and 7 that the use of the PID 
corrector has a great influence on the convergence of the capacity whatever its initial value. This new method 
(EKF combined with PID corrector) characterized by a maximum reliability and minimum response time. 
Regarding the resistance, the results obtained with the new algorithm are the same of those obtained with 
EKF without PID corrector (presented in Figure 8). In the case of resistance, the correction step integrated in 
EKF is enough (Figure 7). The PID corrector proposed in this paper act only on capacitance. 


8. CONCLUSION 

In this paper, a new online identification method applied to supercapacitor is proposed and tested on 
two different supercapacitors. The proposed method is based on extend Kalman filter combined with a PID 
corrector. The supercapacitor was modeled using the equivalent electrical circuit of Zubieta. The proposed 
method allows the online identification of SC capacitance and internal resistance, which are the main 
indicators of SC state-of-health. Unlike other identification method, such as electrochemical impedance 
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spectroscopy or Maxwell test, the proposed method doesn’t require the interruption of system operation. The 
proposed method only requires current and bias voltage measurement, which can be performed online and 
with only two sensors. It has been found that, when the initial capacitance value is far from the real value, the 
extended Kalman filter is unable to accurately determine the correct value of capacitance. For this reason, a 
complementary PID corrector was integrated to EKF in this paper. The capacitance value obtained using the 
new algorithm is very closed to the capacitance obtained using Maxwell test, even when the initial value of 
capacitance is far from the correct value. The absolute errors found are 1.02% and 2.87% for BCAP1500 and 
fBCAP350. Regarding the resistance, the values obtained by the proposed method are also very closed to 
those obtained with EKF. The absolute errors found are 3,71% and 4.02% respectively for BCAP1500 and 
BCAP350. In the case of resistance, the correction step integrated in EKF is enough. In future contributions, 
the identification method proposed in this paper will be applied to other equivalent electrical circuits such as 
multi pore model and simple pore model. The obtained results will be compared to select the most accurate 
model. Then, SC ageing phenomena will be investigated in the goal to evaluate EDLC state of health and 
predict the remaining lifetime using the estimated resistance and capacitance of the SC. 


ACKNOWLEDGEMENTS 
The authors would like to thank Mr Hamid Gualous, the director of LUSAC laboratory. This re- 
search work is sponsored by LUSAC laboratory of the Normandy University, France. 


REFERENCES 

[1] J. R. Miller, and P. Simon, “Materials science: Electrochemical capacitors for energy management,” Science, vol. 
321, no. 5889, pp. 651-652, 2008, doi: 10.1126/science.1158736. 

[2] A. Neffati, M. Guemri, S. Caux, and M. Fadel, “Energy management strategies for multi source systems,” Electric 
Power Systems Research, vol. 102, pp. 42-49, 2013, doi: 10.1016/j.epsr.2013.03.008. 

[3] A. Burke, “Ultracapacitors: why, how, and where is the technology,” Journal of Power Sources, vol. 91, no. 1, pp. 
37-50, 2000, doi: 10.1016/S0378-7753(00)00485-7. 

[4] Q. Zhang, and G. Li, “Experimental study on a semi-active battery-supercapacitor hybrid energy storage system for 
electric vehicle application,’ IEEE Transactions on Power Electronics, vol. 35, no. 1, pp. 1014-1021, Jan. 2020, 
doi: 10.1109/TPEL.2019.2912425. 

[5] B. Yang et al., “Applications of battery/supercapacitor hybrid energy storage systems for electric vehicles using 
perturbation observer based robust control,” Journal of Power Sources, vol. 448, p. 227444, 2020, doi: 
10.1016/).jpowsour.2019.227444. 

[6] Poonam, K. Sharma, A. Arora, and S. K. Tripathi, “Review of supercapacitors: Materials and devices,” Journal of 
Energy Storage, vol. 21, pp. 801-825, 2019, doi: 10.1016/j.est.2019.01.010. 

[7] G. Ren, J. Liu, J. Wan, Y. Guo, and D. Yu, “Overview of wind power intermittency: Impacts, measurements, and 
mitigation solutions,” Applied Energy, vol. 204, pp. 47-65, 2017, doi: 10.1016/j.apenergy.2017.06.098. 

[8] G. Subhashini, R. Abdulla, and T. R. R. Mohan, “Wind turbine mounted on a motorcycle for portable charger,” 
International Journal of Power Electronics and Drive Systems (IJPEDS), vol. 9, no. 4, pp. 1814-1822, 2018, doi: 
ipeds.v9.14.pp1814-1822. 

[9] H. Oussama, O. Abdelkhalek, M. A. Hartani, C. Abdeselem, and M. A. Soumeur, “Proposed energy management 
for a decentralized DC-microgrid based PV-WT-HESS for an isolated community,” International Journal of Power 
Electronics and Drive Systems (IJPEDS), vol. 11, no. 4, pp. 2073, 2020, doi: 10.1159 1/jpeds.v11.14.pp2073-2082. 

[10] M. S. Masaki, L. Zhang, X. Xia, “A hierarchical predictive control for supercapacitor-retrofitted grid-connected 
hybrid renewable systems,” Applied Energy, vol. 242, pp. 393-402, 2019, doi: 10.1016/j.apenergy.2019.03.049. 

[11] R. K6tz and M. Carlen, “Principles and applications of electrochemical capacitors,” Electrochimica Acta, vol. 45, 
no. 15-16, pp. 2483-2498, 2000, doi: 10.1016/S0013-4686(00)00354-6. 

[12] M. H. A. Malek, F. Mustafa, A. M. M. Asry, “A battery-less power supply using supercapacitor as energy storage 
powered by solar,” International Journal of Power Electronics and Drive Systems (IJPEDS), vol. 10, no. 1, pp. 
568, 2019, doi: 10.1159 1/jpeds.v10.11.pp568-574. 

[13] R. German, P. Venet, A. Sari, O. Briat, and J.-M. Vinassa, “Comparison of EDLC impedance models used for 
ageing monitoring,” in 20/2 First International Conference on Renewable Energies and Vehicular Technology, 
Hammamet, Tunisia, pp. 224-229, 2012, doi: 10.1109/REVET.2012.6195275. 

[14] D. B. Dubal, O. Ayyad, V. Ruiz, and P. G. Romero, “Hybrid energy storage: the merging of battery and 
supercapacitor chemistries,’ Chemical Society Reviews, vol. 44, no. 7, pp. 1777-1790, 2015, doi: 
10.1039/c4cs00266k. 

[15] A. Burke, Z. Liu and H. Zhao, “Present and future applications of supercapacitors in electric and hybrid vehicles,” 
2014 IEEE International Electric Vehicle Conference (IEVC), 2014, pp. 1-8, doi: 10.1 JOO/IEVC.2014.7056094. 

[16] X. Luo, J. V. Barreras, C. Chambon, B. Wu, and E. Batzelis, “Hybridizing Lead—Acid Batteries with 
Supercapacitors: A Methodology,” Energies, vol. 14, no. 2, pp. 507, 2021, doi: 10.3390/en14020507. 


Int J Pow Elec & Dri Syst, Vol. 12, No. 3, September 2021 : 1521 — 1534 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 1533 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


[36] 


[37] 
[38] 


[39] 


Z. Song, H. Hofmann, J. Li, J. Hou, X. Han, and M. Ouyang, “Energy management strategies comparison for 
electric vehicles with hybrid energy storage system,” Applied Energy, vol. 134, pp. 321-331, 2014, doi: 
10.1016/j.apenergy.2014.08.035. 

A. Muzaffar, M. B. Ahmed, K. Deshmukh, and J. Thirumalai, “A review on recent advances in hybrid 
supercapacitors: Design, fabrication and applications,” Renewable and Sustainable Energy Reviews, vol. 101, pp. 
123-145, 2019, doi: 10.1016/j.rser.2018.10.026. 

B. Sami, N. SIhem, S. Gherairi, and C. Adnane, “A Multi-Agent System for Smart Energy Management Devoted to 
Vehicle Applications: Realistic Dynamic Hybrid Electric System Using Hydrogen as a Fuel,” Energies, vol. 12, no. 
3, pp. 474, 2019, doi: 10.3390/en12030474. 

R. Kotz, P. W. Ruch, and D. Cericola, “Aging and failure mode of electrochemical double layer capacitors during 
accelerated constant load tests,” Journal of Power Sources, vol. 195, no. 3, pp. 923-928, 2010, dot: 
10.1016/).jpowsour.2009.08.045. 

H. Gualous, R. Gallay, M. Al Sakka, A. Oukour, B. T.-Ighil, and B. BOudart, “Calendar and cycling ageing of 
activated carbon supercapacitor for automotive application,” Microelectronics Reliability, vol. 52, no. 9-10, pp. 
2477-2481, 2012, doi: 10.1016/j.microrel.2012.06.099. 

A. Eddahech, O. Briat, N. Bertrand, J. Y. Delerage, and J.-M. Vinassa, “Behavior and state-of-health monitoring of 
Li-ion batteries using impedance spectroscopy and recurrent neural networks,” International Journal of Electrical 
Power & Energy Systems, vol. 42, no. 1, pp. 487-494, 2012, doi: 10.1016/j.14jepes.2012.04.050. 

F. Faroog, A. Khan, S. J. lee, M. M. Nadeem, and W. Choi, “A multi-channel fast impedance spectroscopy 
instrument developed for quality assurance of super-capacitors,” Energies, vol. 14, no. 4, pp. 1139, 2021, dot: 
10.3390/en 14041139. 

A. Nadeau, G. Sharma and T. Soyata, “State-of-charge estimation for supercapacitors: A Kalman filtering 
formulation,” in 20/4 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), 2014, 
pp. 2194-2198, doi: 10.1109/ICASSP.2014.6853988. 

C. J. Chiang, J. L. Yang, and W. C. Cheng, “Temperature and state-of-charge estimation in ultracapacitors based on 
extended Kalman filter,’ Journal of Power Sources, vol. 234, pp. 234-243, 2013, doi: 
10.1016/j.jpowsour.2013.01.173. 

Z. Shi, F. Auger, E. Schaeffer, P. Guillemet and L. Loron, “Interconnected observers for online supercapacitor 
ageing monitoring,” in JECON 2013 - 39th Annual Conference of the IEEE Industrial Electronics Society, 2013, 
pp. 6746-6751, doi: 10.1109/IECON.2013.6700249. 

A. El Mejdoubi, H. Chaoui, H. Gualous, and J. Sabor, “Online Parameter Identification for Supercapacitor State-of- 
Health Diagnosis for Vehicular Applications,” IEEE Transactions on Power Electronics, vol. 32, no. 12, pp. 9355- 
9363, Dec. 2017, doi: 10.1109/TPEL.2017.2655578. 

H. Chaoui, A. El Mejdoubi, A. Oukaour, and H. Gualous, “online system identification for lifetime diagnostic of 
supercapacitors with guaranteed stability,” IEEE Transactions on Control Systems Technology, vol. 24, no. 6, pp. 
2094-2102, Nov. 2016, doi: 10.1109/TCST.2016.2520911. 

P. Saha, and M. Khanra, “Online estimation of state-of-charge, state-of-health and temperature of supercapacitor,” 
in 2020 IEEE International Symposium on Circuits and Systems (ISCAS), 2020, pp. 1-5, dot: 
10.1 109/ISCAS45731.2020.9180494. 

Z. Peng, and L. Jikai, “On new UAV flight control system based on Kalman & PID,” in 2011 2nd International 
Conference on Intelligent Control and Information Processing, 2011, pp. 819-823, doi: 
10.1 1O9/ICICIP.201 1.6008362. 

S. Wakitani, H. Nakanishi, Y. Ashida, and T. Yamamoto, “Study on a Kalman Filter based PID Controller,” IFAC- 
PapersOnLine, vol. 51, no. 4, pp. 422-425, 2018, doi: 10.1016/j.ifacol.2018.06.131. 

B. H. Prasetio, “Ensemble Kalman filter and PID controller implementation on self balancing robot,” in 20/5 
International Electronics Symposium (IES), 2015, pp. 105-109, doi: 10.1109/ELECS Y M.2015.7380823. 

F. Naseri, E. Farjah, T. Ghanbari, Z. Kazemi, E. Schaltz and J. Schanen, “Online parameter estimation for 
supercapacitor state-of-energy and state-of-health determination in vehicular applications,” JEEE Transactions on 
Industrial Electronics, vol. 67, no. 9, pp. 7963-7972, Sept. 2020, doi: 10.1109/TIE.2019.2941151. 

A. Soualhi et al., “Heath monitoring of capacitors and supercapacitors using the neo-fuzzy neural approach,” IEEE 
Transactions on Industrial Informatics, vol. 14, no. 1, pp. 24-34, Jan. 2018, doi: 10.1109/TH.2017.2701823. 

S. Zhao, F. Blaabjerg, and H. Wang, “An Overview of Artificial Intelligence Applications for Power Electronics,” 
IEEE Transactions on Power Electronics, vol. 36, no. 4, pp. 4633-4658, April 2021, doi: 
10.1109/TPEL.2020.3024914. 

Y. Diab, P. Venet, H. Gualous, and G. Rojat, “Self-discharge characterization and modeling of electrochemical 
capacitor used for power electronics applications,” IEEE Transactions on Power Electronics, vol. 24, no. 2, pp. 
510-517, Feb. 2009, doi: 10.1109/TPEL.2008.20071 16. 

Z. Bououchma, J. Sabor, and H. Aitbouh, “New electrical model of supercapacitors for electric hybrid vehicle 
applications,” in Materials Today: Proceedings, vol. 13, pp. 688-697, 2019, doi: 10.1016/).matpr.2019.04.029. 

L. Zubieta, and R. Bonert, “Characterization of double-layer capacitors for power electronics applications,” IEEE 
Transactions on Industry Applications, vol. 36, no. 1, pp. 199-205, Jan.-Feb. 2000, doi: 10.1109/28.821816. 

L. Zhang, Z. Wang, X. Hu, F. Sun, and D. G. Dorrell, “A comparative study of equivalent circuit models of 
ultracapacitors for electric vehicles,” Journal of Power Sources, vol. 274, pp. 899-906, 2015, doi: 
10.1016/j.jpowsour.2014.10.170. 


Online diagnosis of supercapacitors using extended Kalman filter ... (Zoubida Bououchma) 


1534 O ISSN: 2088-8694 


[40] R. German, A. Hammar, R. Lallemand, A. Sari, and P. Venet, “Novel experimental identification method for a 
supercapacitor multipore model in order to monitor the state of health,” IEEE Transactions on Power Electronics, 
vol. 31, no. 1, pp. 548-559, Jan. 2016, doi: 10.1109/TPEL.2015.2408457. 

[41] G. Krishnan, S. Das, and V. Agarwal, “An online identification algorithm to determine the parameters of the 
fractional-order model of a supercapacitor,” IEEE Transactions on Industry Applications, vol. 56, no. 1, pp. 763- 
770, Jan.-Feb. 2020, doi: 10.1109/TIA.2019.2947235. 

[42] L. Zhang, X. Hu, Z. Wang, F. Sun, and D. G. Dorrell, “A review of supercapacitor modeling, estimation, and 
applications: A control/management perspective,” Renewable and Sustainable Energy Reviews, vol. 81, pp. 1868- 
1878, 2018, doi: 10.1016/j.rser.2017.05.283. 

[43] A. Oukaour, B. Tala-Ighil, M. AlSakka, H. Gualous, and R. Gallay, “Calendar ageing and health diagnosis of 
supercapacitor,” Electric Power Systems Research, vol. 95, pp. 330-338, 2013, doi: 10.1016/j.epsr.2012.09.005. 

[44] R. Toscano, and P. Lyonnet, “A Kalman Optimization Approach for Solving Some Industrial Electronics 
Problems,” IEEE Transactions on Industrial Electronics, vol. 59, no. 11, pp. 4456-4464, Nov. 2012, dot: 
10.1109/TIE.2011.2169637. 

[45] B. Gou, “Jacobian matrix-based observability analysis for state estimation,” IEEE Transactions on Power Systems, 
vol. 21, no. 1, pp. 348-356, Feb. 2006, doi: 10.1109/TPWRS.2005.860934. 


BIOGRAPHIES OF AUTHORS 


RA 


~e 


aw 


“a. d if 


Zoubida Bououchma was born in Morocco in 1992. She received the Engineer degree in 
Electromechanical Engineering in 2015 from the Ecole Nationale Supérieure d’Arts & Métiers 
of Meknes (ENSAM-Meknes), Morocco. She is currently a Ph.D. student in electrical 
engineering in ENSAM-Meknes. His research interests are in the area of energy storage, more 
specifically, supercapacitors and Lithium-ion batteries, State-of-Health and State-of-Charge 
diagnosis and aging estimation. 








Jalal Sabor received the PhD degree in engineering science from the Institut National des 
Sciences Appliquées (INSA), Rouen, France in 1995. He is currently a full Professor of 
industrial computer science at the Ecole Nationale Supérieure d’Arts & Métiers (ENSAM) 
Moulay Ismail university, Meknes, Morocco. He is a member of the LSMI Laboratory, He is 
also the research team control steering and supervision systems head. His actual main research 
interests concern intelligent management of energy, smart grid, control and supervision systems, 
Architecture Based on Multi Agents Systems and fuzzy logic. 





Int J Pow Elec & Dri Syst, Vol. 12, No. 3, September 2021 : 1521 — 1534 


